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Abstract 

This is a follow-up to our earlier work on the energies and radial distributions of heavy-light 
mesons. The heavy quark is taken to be static (infinitely heavy) and the light quark has a mass 
about that of the strange quark. We now concentrate on the energies of the excited states with 
higher angular momentum and with a radial node. A new improvement is the use of hypercubic 
blocking in the time direction. 

The calculation is carried out with dynamical fermions on a 16 3 x 32 lattice with a lattice spacing 
a as 0.1 fm generated using a non-perturbatively improved clover action. 

In nature the closest equivalent of this heavy-light system is the B s meson, which allows us to 
compare our lattice calculations to experimental results (where available) or to give a prediction 
where the excited states, particularly P-wave states, should lie. We pay special attention to the 
spin-orbit splitting, to see which one of the states (for a given angular momentum L) has the lower 
energy. An attempt is made to understand these results in terms of the Dirac equation. 
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I. MOTIVATION 



There are several advantages in studying a heavy-light system on a lattice. First of all, 
the lattice calculations are relatively easy to do, and it allows us to do these with QCD 
from first principles. Our meson is much more simple than in true QCD: one of the quarks 
is static (infinitely heavy) with the light quark "orbiting" it. This makes it very beneficial 
for modelling. On the lattice an abundance of data can be produced, and we know which 
state we are measuring. In contrast, the physical states can be a mixture of two or more 
configurations, but on the lattice this complication is, mostly, avoided. Even so, our results 
on the heavy-light system can still be compared to the B s meson experimental results. 



II. MEASUREMENTS AND LATTICE PARAMETERS 

We measure the energies of both angular and first radial excitations of heavy-light mesons. 
Since the heavy quark spin decouples from any description of the configurations we may label 
the states as Ly- = L ± ~, where L is the orbital angular momentum and ±| refers to the 
spin of the light quark. 

The measurements are done on 16 3 x 32 lattices using two degenerate quark flavours. 
The lattice configurations were generated by the UKQCD Collaboration using lattice action 
parameters (3 = 5.2, csw = 2.0171 and three different values for the hopping parameter k 
(see Table HJ. The three different lattices are referred to here as "DF3", "DF4" and "DF5". 
Each of them has a slightly different lattice spacing (a) and a different light quark mass 
(m q ). Our main results are measured on the "DF3" lattice, because the light quark mass is 
very close to the strange quark mass. More details of the lattice configurations used in this 
study can be found in Refs. . Because our light quarks are heavier than true u and d 
quarks, we have ranging from 730 MeV ("DF3") to 400 MeV ("DF5"). Two different 
levels of fuzzing (2 and 8 iterations of conventional fuzzing) are used in the spatial directions 
to permit a cleaner extraction of the excited states. 
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III. 2-POINT CORRELATION FUNCTION 

The 2-point correlation function (see Fig. [T]) is denned as 

C 2 (T) = (P t TG q (x,t + T,t)P t+T ^U Q (x,t,t + T)) , (1) 

where U®(x,t,t + T) is the heavy (infinite mass)-quark propagator and G q (x,t + T,t) the 
light anti-quark propagator. P t is a linear combination of products of gauge links at time t 
along paths P and T defines the spin structure of the operator. The (...) means the average 
over the whole lattice. A detailed discussion of lattice operators for orbitally excited mesons 
can be found in [31]. In this study, the same operators are used as in |4j. The energies (mj) 
and amplitudes (cij) are extracted by fitting the C2 with a sum of exponentials, 

N m ax 

[C2{T)] h , h « Yl a U^~ m ' T ^f2, where A max = 2 - 4, T < 14. (2) 
i=i 

The fit is a simple least squares fit. In most of the cases 3 exponentials are used to try 
to ensure the first radially excited states are not polluted by higher states. Also 2 and 4 
exponential fits are used to cross-check the results wherever possible. Indices fi and f 2 
denote the amount of fuzzing used at the vertices and both of them take two values, fx =F1, 
F2 and f'2 =F1, F2, where (Fl=2 iterations and F2=2+6 iterations). For S and P_ states 
we have alternative operators (see j^j]), so we get a 5 by 5 matrix (5 paths, because one 
operator has two choices, Fl and F2, and the other operator has three choices, local, Fl and 
F2) instead of just a 2 by 2 matrix (2 paths) given by the fuzzing choices. 





K 


r /a 


a [fm] (approx.) 


mq/m s (approx.) 


r m n 


No. of configs. 


DF3 
DF4 
DF5 


0.1350 
0.1355 
0.1358 


4.754(40)^ 
5.041(40);? 
5.32(5) 


0.110 
0.104 
0.099 


1.1 

0.6 
0.3 


1.93(3) 
1.48(3) 
1.06(3) 


160 
119 
139 



TABLE I: Lattice parameters (from [l|). The Sommer scale parameter can be taken to be ro = 
0.525(25) fm and m s is the s quark mass. 
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IV. SMEARED HEAVY QUARK 



We introduce two types of smearing in the time direction to get a better noise to signal 
ratio. The first type is APE type smearing, where the original links in the time direction 
are replaced by a sum over the six staples that extend one lattice spacing in the spatial 
directions (in Fig. [2] on the left). This smearing is called here "sum6" for short. We use the 
notation "plain" to refer to the original Eichten-Hill point static source construction. 

To smear the heavy quark even more we then use hypercubic blocking (first introduced 
by Hasenfratz and Knechtli in [5|), again only for the links in the time direction (in Fig. [2] 
on the right). Now the staples (the red dashed lines in Fig. [2]) are not constructed from the 
original, single links, but from staples (the blue dash-dotted lines in Fig. [2]). In more detail, 
we first construct the links 



V, 



I, p, ; v p 



Proj 



SU(3) 



;i-a 3 )^ + y J2 U ^ U -< 

±t] + p,v,p 



i+V,V u i+p,ri 



(3) 



where Ui jfl is the original thin link at location i and direction /i. Note that there are no 
staples in directions v or p. We then construct "fat" links 



Vi 



I, p ; v 



Proj 



SU(3) 



q>2 \ ^ — — — + 



(4) 



where index v indicates that the link is not decorated with staples in that direction. The 
last step is 



Vi 



Proj 



SU(3) 



ai)U t 



i, p. 



Q'l 
If 



Vi,v\pyi+v,p,-,vVi-\ 

±v + p 



(5) 




FIG. 1: Two-point correlation function. 
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FIG. 2: (Color online) APE smearing in the time direction (on the left) and hypercubic blocking 
(on the right). 

where the "fat" links are again used to construct the new links. The values c*3 = 0.5, «2 = 1 
and oil = 1 are used in this study, because this choice was found to be very good in reducing 
the noise to signal ratio in {(J. Note that = 1.0, «2 — and a± = would give the 
"sum6" smearing. Hypercubic blocking takes into account the links within a "hypercube" 
(the edges of the "cube" are 2a in spatial directions but only one lattice spacing in the 
time direction). This smearing is called here "hyp" for short. The "plain" configurations 
do not have smearing in the time direction. Smearing the heavy quark improves the noise 
to signal ratio, which can be seen in Figs. [3H3 The figures show the standard deviation to 
signal ratio for the largest component of C2, which is F2F2, for 160 lattice configurations 
(lattice "DF3"). In all cases the "plain" signal is clearly inferior to the "sum6" and "hyp" 
signals, whereas the "hyp" signal is also better than the "sum6" signal. This latter difference 
wouldbe more apparent in a non- logarithmic scale. Lattices "DF4" and "DF5" show similar 
trends as the "DF3" lattice. 
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FIG. 3: (Color online) Noise (standard deviation) to signal ratio: S-wave 2-point correlation func- 
tion C*2 for the lattice "DF3" . Note the logarithmic scale. 
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FIG. 4: (Color online) Noise (standard deviation) to signal ratio as in Fig. [3] but for the P_ case. 



6 



a 

CO 



10 



1 r 



~ 0.1 =r 



g o.oib- 

'3 

z 



0.001 r 



0.0001, 



O DF3plain 

□ DF3sum6 
DF3hyp 



9 



9 9 



@ 8 



o o 



o n 5 o 
^ o □ 

O 50 



^~5" 
□ 

o 1 







10 



12 



14 



16 



FIG. 5: (Color online) Noise (standard deviation) to signal ratio for the P+ case. 
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FIG. 6: (Color online) Noise (standard deviation) to signal ratio for the D_ case. 
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FIG. 7: (Color online) Noise (standard deviation) to signal ratio for the D + case. 



L± 


mi 


m 2 


m 3 


«1,F1 


Ol,F2 


02,F1 


G2,F2 


03,F1 


03,F2 


x 2 

d.o.f. 


S 


0.527(5) 


0.98(2) 


1.40(2) 


0.297(11) 


0.44(2) 


1.24(5) 


1.56(4) 


1.81(5) 0.03(11) 


21/24 


P 


0.766(14) 


1.29(3) 


1.52(3) 


0.65(5) 


0.93(7) 


3.4(3) 


2.2(2) 


-1.5(6) 


3.4(6) 


19/24 


P+ 


0.76(2) 


1.28(3) 


1.46(2) 


0.62(7) 


0.88(10) 


3.6(4) 


3.14(15) 


-2.8(5) 


2.6(5) 


11/24 


D+- 


1.10(8) 


1.46(5) 


1.66(6)* 


1.0(3) 


2.0(5) 


5.3(8) 


1.8(9) 


-0.6(17) 


5.1(13) 


22/24 


D 


1.01(3) 


1.52(2) 


1.67(4)* 


1.10(14) 


2.0(2) 


7.2(5) 


1.5(12) 


-1.4(15) 


7.4(8) 


36/27 


D+ 


1.06(2) 


1.558(8) 


1.80(2)* 


1.14(9) 


2.09(11) 


6.52(6) 


-0.5(5)0 


1.3(7) 


7.6(2) 


39/27 


F+- 


1.20(2) 


1.658(5) 


1.892(13)* 


0.77(8) 


1.88(10) 5.30(13) 


-1.1(5) 


1.3(7) 


6.0(3) 


37/27 



TABLE II: Two-point correlation function fits (equation |2J) for "DF3hyp". 2 path fit results 
are shown for all states to make comparisons easier, even though our best fits for the S and 
P_ states are 5 path fits. In some cases (entries marked with an asterisk) Bayesian ideas are 
used (see section IV B[) . The errors on the parameters were obtained by bootstrapping the lattice 

2 

configurations and repeating the fit 100 times. The mj are in lattice units. Note that the is 
larger for the D-wave and F-wave states than for the S- and P-wave states. 
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FIG. 8: (Color online) Effective mass plot for the S-wave F2F2 correlations. The line labelled "fit" 
(here and in the other effective mass plots) shows the lowest energy obtained from the fit in Eq. [2] 
for the lattice "DF3hyp" . Only 2 paths are used in the fit shown here to make comparisons easier, 
although the best fit for the S-wave state uses all 5 paths. The thickness of the line indicates the 
error. As expected from a variational argument, this fit to all data naturally gives a somewhat 
smaller mass than a fit to F2F2 alone. 

V. ENERGY SPECTRUM 

The energies are obtained from the fit in equation [2] — see Table HT1 for the results for the 
lattice "DF3hyp" . The m« are in lattice units. However, due to the presence of an unknown 
(but L± and % independent) self energy in each m,, only the differences m*(L±) — mi(S) are 
relevant. The ground state energy from the 2-path fit for a given state is compared with the 
effective mass in Figs. I81TT21 To illustrate how plateaux develop with T, In [C%(T) / 'C^iT + 1)] 
is shown for the largest component, F2F2, for the same lattice "DF3hyp". It is seen that 
for all states we get a plateau that agrees nicely with the fit result that uses all components, 
F1F1, F1F2 and F2F2. The errors are large for large T values, and the data points for the 
highest T values are not shown if the errors render them insignificant. In the fits this is 
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FIG. 9: (Color online) Effective mass plot for the P_ F2F2 correlations. Only 2 paths are used in 
the fit shown here to make comparisons easier, although the best fit for this state uses all 5 paths. 
Other details as in Fig. El 

under better control, because we fit C^T) and not the ratio, and all fuzzing combinations 
are used (i.e. more data are used). We can thus use data up to T = 15 in the fits. The 
fit shown in these figures is only to the 2 path data, to make comparisons easier. When 
extracting the energy of a state all 5 paths are used for S and P_. For other states only 2 
path data are available. 

The resulting energy spectra from the fit (Eq. [5]) for different lattices are shown in 
Figs. [T^rfTol — see also Tables [TTTHVl With the lattice "DF3", for most states, using different 
smearing for the heavy quark does not seem to change significantly the energy differences 
with respect to the IS energy — the exceptions being the P + and excited D + _ states. 
Different smearings should only give the same results in the continuum limit, so it is under- 
standable that at a fixed lattice spacing the results may differ. Unfortunately, all our lattices 
have approximately the same lattice spacing (about a = 0.1 fm, see Table H]), and we can not 
go to the continuum limit properly. However, we can use the results from different smearings 
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FIG. 10: (Color online) Effective mass plot for the P+ F2F2 correlations. Other details as in Fig.0 
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FIG. 11: (Color online) Effective mass plot for the D_ F2F2 correlations. Other details as in Fig. El 
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FIG. 12: (Color online) Effective mass plot for the D + F2F2 correlations. Other details as in Fig. 



to give a rough estimate of the systematic error. Because studying the noise to signal ratio 
(Figs. EHZ]) shows that the "plain" configurations are clearly inferior to the configurations 
that are smeared in the time direction, we use "hyp" and "sum6" configurations to get our 
main results. The reason for quoting energies in units of tq is to avoid the 5% uncertainty 
in ro = 0.525(25) fm. The uncertainties in r^/a are much smaller. To get energies in GeV 
then requires an additional factor of 0.38(2). 

The energy of the D + _ state had been expected to be near the spin average of the D_ 
and D + energies, but it turns out to be a poor estimate of this average. Therefore, it 
is not clear to what extent the F + _ energy is near the spin average of the two F-wave 
states, as was originally hoped. The F- and D-wave results for different smearings for the 
lattices "DF4" and "DF5" are somewhat more scattered than for the "DF3" lattice (i.e. 
the systematic errors are larger), but otherwise the same features are seen for all three light 
quark masses. We are most interested in the "DF3" lattice results, because that is closest 
to the B s meson (the light quark mass on this lattice being close to the s quark mass). One 
interesting observation is, that the energy spectrum is close to being dependent on L alone. 
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FIG. 13: (Color online) Effective mass plot for the F_| F2F2 correlations. Other details as in 

Fig. I 

For example our preferred configurations, "DF3hyp", show an approximate linear rise in 
excitation energy with L (up to F-wave) as ~0.4L GeV. A similar linear rise is usually seen 
in Regge or string models. In contrast, the 2S state is seen to be almost degenerate with 
the ID states, as in a simple harmonic oscillator. A L(L+1) term can be added to the linear 
ansatz to get a better fit — more precisely, 0.34L+0.04L(L+1) gives a good overall fit to 
the four energy differences up to D-waves. 

Our earlier results (for the "plain" configurations used in this study as well as for some 
other unquenched and quenched configurations) can be found in Refs. . Because dy- 
namical fermions are used in this study, we can not be absolutely certain that the lattice 
states are pure quark-anti-quark states. However, our radial distribution measurements 8] 
support the assumption that the states are ordinary meson states: the radial distributions 
of the lowest lying states are not broad, as would be the case if the states were molecules, 
and the first radial excitations of S- and P-wave states have one node at short distances 
(approximately at 0.30-0.35 fm). 
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FIG. 14: (Color online) Energy spectrum of the heavy-light meson using lattice "DF3" in units of 

ro . Here L+(_) means that the light quark spin couples to angular momentum L giving the total 

j = L± 1/2. The 2S is the first radially excited L= state. The D_| is a mixture of the D_ and 

D + states, and likewise for the Fj Energies are given with respect to the S-wave ground state 

n 

(IS). Here r /a = 4.754(40)1^ (from [1]). The error bars shown here contain only the statistical 
errors on the lattice energy fits. 
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FIG. 15: (Color online) Energy spectrum of the heavy-light meson using lattice "DF4". See 
Fig. [H for details. Here r /a = 5.041(40)+ ( j ) (from [1]). The error bars shown here contain only 
the statistical errors on the lattice energy fits. 
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FIG. 16: (Color online) Energy spectrum of the heavy-light meson using lattice "DF5". See Fig.[T4] 
for details. Here r^/a = 5.32(5) (from |l[). The error bars shown here contain only the statistical 
errors from the lattice energy fits. 





DF3plain 


DF3sum6 


DF3hyp 


nL± 


DF3plain 


DF3sum6 


DF3hyp 


1S 


3.55(6) 


2.724(14) 


2.520(10) 


2S 


5.1(2) 


4.64(10) 


4.59(8) 


1P_ 


4.74(3) 


3.79(4) 


3.62(3) 


2P_ 


7.13(5) 


5.98(10) 


5.97(10) 


1P+ 


5.15(11) 


4.12(8) 


3.63(10) 


2P+ 


7.80(9) 


6.62(8) 


6.1(2) 


lD+_ 




5.48(9) 


5.25(14) 


2D+_ 




6.88(4) 


7.05(7) 


1D_ 


6.1(2) 


5.10(7) 


4.79(13) 


2D_ 


8.44(14) 


7.50(4) 


7.23(11) 


1D+ 


6.2(2) 


5.23(9) 


5.06(7) 


2D + 


8.38(15) 


7.59(5) 


7.41(4) 


1F+^ 


7.21(12) 


6.10(11) 


5.69(8) 


2F+_ 


9.16(3) 


8.18(3) 


7.88(2) 



TABLE III: Heavy-light meson energies on the lattice in units of tq for "DF3". The uncertainty 
due to the statistical error on r^/a [rp/o = 4.754(40)lg , from [1]] is small (less than 1% ), and is 
not taken into account in the error estimates. The n denotes the radial excitation and n— 1 gives 
the number of nodes in the wavefunction of the state. The dash means that no reliable fit can 
be found. The results on different lattices can not be compared directly (only energy differences 
can be compared) due to different self energies. The "DF3hyp" results are the same as mi, m,2 in 
Table HI] but expressed in different units. Also the IS, 1P_ are now from 5 path fits. 
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2.45(2) 
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4.77(11) 
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3.43(5) 
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6.03(12) 


5.91(11) 
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3.37(16) 
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4.4(3) 


2D + 


9.16(8) 


7.9(2) 


7.1(2) 


lF+_ 


6.8(4) 


6.24(10) 


5.7(2) 


2F + _ 


9.51(11) 


8.62(2) 


8.33(10) 



TABLE IV: Heavy-light meson energies on the lattice in units of ro for "DF4" . The uncertainty 
due to the statistical error on ro/a [ro/a = 5.041 (40)1^0 > f rom is small (less than 1% ), and is 
not taken into account in the error estimates. Other comments as in Table IIIII 
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2.70(3) 


2.46(3) 


2S 


5.70(6) 
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4.63(14) 
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4.75(4) 
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6.11(10) 
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4.01(15) 
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6.3(2) 


1D+- 


6.6(3) 






2D+_ 


8.71(13) 






1D_ 


6.86(12) 


4.7(3) 


4.6(2) 


2D_ 


9.51(5) 


7.5(2) 


7.4(2) 


1D+ 


7.0(2) 


5.2(3) 


4.7(5) 


2D + 


9.61(8) 


8.2(2) 


7.6(3) 


lF + _ 


7.7(2) 


6.72(10) 


6.3(2) 


2F + _ 


10.05(6) 


9.01(3) 


8.6(3) 



TABLE V: Heavy-light meson energies on the lattice in units of ro for "DF5" . The uncertainty 
due to the statistical error on r$/a [ro/a = 5.32(5), from [l]] is small (less than 1% ), and is not 
taken into account in the error estimates. Other comments as in Table IIIII 

To check how the results depend on the light quark mass we plot the energies (i.e. energy 
differences with respect to the IS state) as a function of the pion mass squared, (r m n ) 2 . As 
can be seen in Figs. [T7H2DI for P- and D-wave states the dependence on the light quark mass 
is not strong. We also compare our results to other static-light meson lattice calculations in 
Fig. EU It is seen that the P-wave results do not change much between the different lattices, 
but the 2S-1S and D-wave energy differences vary a lot. However, since the lattices, quark 
masses and lattice spacings are different, the results should only agree in the continuum 
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FIG. 17: (Color online) The energy difference ro[E(lP_)— E(1S)] as a function of the pion mass 
squared. There is a slight dependence on the light quark mass. Here and in the following figures 
for P + , D and D + the results, from left to right, are from lattices "DF5", "DF4" and "DF3" , 
respectively. 



limit. A set of similar lattices with different lattice spacings would be needed to check this. 



A. Interpolation to the b quark mass 

Even though we can not go to the continuum limit, we feel that it is worth-while to try to 
predict where the B s meson excited states lie. To obtain the predictions of the excited state 
energies, we can now interpolate in l/rng, where ijiq is the heavy quark mass, between the 
"DF3" heavy-light lattice calculations and D s meson experimental results, i.e. interpolate 
between the static quark (tuq = oo) and the charm quark (juq = m c ). Here we, of course, 
have to assume that the measured D s meson states are simple quark-anti-quark states. This 
is not necessarily true: for example the mass of the _D* (2317) is much lower than what is 
predicted by conventional potential models, and it has thus been proposed that it could be 
either a four quark state, a DK molecule or a Dn atom. A short review on meson excited 
state spectroscopy and the puzzles in interpreting the results is given in Ref. [ljj. However, 
the inclusion of chiral radiative corrections could change the potential model predictions 



17 



0.8 : 
0.6 - 



V sum6 
> hyp 







( r n m ) 

v 7T 



FIG. 18: (Color online) The energy difference ro[E(lP+)— E(1S)] as a function of the pion mass 
squared. The dependence on the light quark mass is weak, whereas there is a manifest difference 
between "sum6" and "hyp" configurations. 
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FIG. 19: (Color online) The energy difference ro[E(lD_)— E(1S)] as a function of the pion mass 
squared. The dependence on the light quark mass is weak, and the differences between the "sum6" 
and "hyp" smearings are small — especially for "DF3" . 
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FIG. 20: (Color online) The energy difference ro[E(lD+)— E(1S)] as a function of the pion mass 
squared. The dependence on the light quark mass is weak, and the differences between the "sum6" 
and "hyp" smearings are small — especially for "DF3" . 

considerably [ijj . In the following we assume that the states are the usual quark-anti-quark 
states. 

We use linear interpolation, i.e. 



AE 



m. 



A + B— + CF r 

m Q m Q 



(6) 



Here A, B and C are fit parameters and Fj = 2[J(J + 1) — j q (j q + 1) — sq(sq + 1)], where 
J is the total angular momentum, sq = 1/2 is the heavy quark spin and j q is the combined 
spin and orbital angular momentum L of the light quark (see Table IVll) . The interpolation 
procedure is shown in Fig. [221 Note that the linear interpolation works perfectly for the 1~ 
S-wave state, where the experimental energies are known for both B s and D s mesons, and 
the lattice result (zero) is simply because the two IS states are automatically degenerate at 
niQ = oo. 

Our predictions of the energy differences m(lP) — m(lS) for the B s meson are given in 
Table IVIII For our preferred lattice "DF3hyp" these agree very well with the experimental 
measurements of the energies of the 1 + and 2 + P-wave states. There we predict that the two 
lowest P-wave states lie a few MeV below the BK and B*K thresholds (minus the IS state 
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FIG. 21: (Color online) Comparison of different static-light lattice results. On the left at a = 
0.0855 fm are the results from European Twisted Mass Collaboration [jj], at a = 0.11 fm our 
results, at a = 0.16 fm BGR Collaboration's results 101 and on the right at a = 0.17 fm TrinLat 

n 

group's results [15|. A set of similar lattices with different lattice spacings is needed for going to 
the continuum limit. 





J P L Fj 


J P L Fj 


0+ -3 
1+ +1 


0+ 1 -3 
1+ 1 +1 


1+ 1 -5 
2+ 1 +3 



TABLE VI: Coefficients Fj = 2[J(J + 1) - j q (j q + 1) - s Q (s Q + 1)] (equation ED . 

energy) at 406 and 452 MeV respectively. We show the "DFsum6" results for comparison. 

As for other excited states, BaBar and Belle observed two new states, _D*j (2860) and 
£^(2700), in 2006 HQ- The J p quantum numbers of the D* sJ (28Q0) can be + , 1", 2+ 
etc., so it could be a radial excitation of the D* (2317) or a J p = 3~ D-wave state. The first 
interpretation is rather popular, but our lattice results favor the D-wave J p = 3 _ assignment 
in agreement with Colangelo, De Fazio and Nicotri 18) . (The slope of the interpolating line 
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FIG. 22: (Color online) Interpolation to the b quark mass. The ratio m c /mb is taken to be 0.30(2) 



(from 11]; shown by the vertical band). The D s meson experimental results are from 



circles), and the B s meson experimental results are from 



I 1 1 (blue circles) and [3| (green triangles ) 



111 (blue 



Our results (using "DF3hyp" configurations) are marked with red squares. 





DF3hyp 


DF3sum6 


experiments 


0+ 


393 ± 9 MeV 


384 ± 10 MeV 




1+ 


440 ± 9 MeV 


432 ± 10 MeV 




1+ 


466 ± 25 MeV 


538 ± 21 MeV 


463 ± 1 MeV 


2+ 


482 ± 25 MeV 


551 ± 21 MeV 


473 ± 1 MeV 



TABLE VII: Our predictions for B s meson mass differences, M(B*)-M(B S ), for the P-wave states. 
The uncertainty in the ratio m c /mb is not taken into account in the error estimates. The experi- 
mental results are from 121. 
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FIG. 23: (Color online) Interpolation to the b quark mass for "DF3hyp" lattice: higher excited 
states. The lines illustrate what the interpolation would look like, if the D s meson states were D- 
wave states. The experimental results are from [lrl [l7]. Interpolating to m c /m;, predicts D-wave 
J p = 1~, 3- at 817(31) and 932(18) MeV respectively. 



would be very steep, if the -D*j(2860) is a radial excitation of the -D* (2317).) Interpolation 
then predicts a D-wave J p = 3" B s state at 932(18) MeV. In addition, the £)* J (2700) could 
be a radially excited S-wave state or a D-wave J p = 1' state. If the latter identification is 
assumed, then a D-wave J p = 1~ B s state at 817(31) MeV is expected (see Fig. [231 . 



B. Bayesian ideas 

In some cases, using 3 exponentials to fit the Gi data does not work very well. In Table HTl 
these cases are marked with an asterisk. Since these fits are not as good or stable as one 
would hope, we introduce some Bayesian ideas and use prior knowledge of the energies to 
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FIG. 24: (Color online) Mass differences A77732 = rn.3 — m<i for the "DF3hyp" data in lattice units. 
This seems to be almost constant for angular momentum L> 1 (given the sizeable errors on the 
data). The lines give the Am^ = 0.234(46) that is used in this study. Looking at the "DF3sum6" 
data gives a very similar picture (not shown here) and an estimate of A77732 = 0.207(54). See 
section IV Bl for details. 



nL± 


777,3 3 exp 


^3, prior 


777,3 Bayes 


ni2 3 exp 


TO2 Bayes 


P+ 


1.46(2) 


1.52(6) 


1.471(15) 


1.28(3) 


1.29(3) 


D_ 


1.63(5) 


1.74(6) 


1.67(4) 


1.51(3) 


1.52(2) 


D+ 


1.84(12) 


1.79(5) 


1.80(2) 


1.559(11) 


1.558(8) 


F+- 


1.96(43) 


1.89(5) 


1.89(1) 


1.66(2) 


1.657(5) 



TABLE VIII: Comparison of the m3 jPr ior with the results from the full 3 exponential fit and the 
Bayesian fit for "DF3hyp" configurations. See Section |VB] for definition of m.3 iPr ior- " m 3 3 exp" and 
"7713 Bayes" are the results of a full 3 exponential fit and a Bayesian (fixed 771.32) fit, respectively. 
Likewise for the 7712- The P+ "Bayes" fit is merely to check that the Bayesian ideas work well and 
does not restrict the analysis too much. 

constrain the fit, or rather to guide the fit in the right direction. The third mass, m 3 , (which 
would be the mass of the second radial excitation, if there was no pollution from higher 
states) is restricted to be in the range m 3) prioi . ± Am 3> prioi . by adding a term 



(m 3 - m 3 , pr ior) ' 

(Am 3) prior) 2 
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(7) 



nL± 


7713 3 exp 


^3, prior 


7773 Bayes 


rrt2 3 exp 


mi Bayes 


P+ 


1.59(2) 


1.61(5) 


1.59(2) 


1.39(2) 


1.393(15) 


D_ 


1.80(9) 


1.79(5) 


1.79(3) 


1.578(11) 


1.577(8) 


D+ 


2.1(2) 


1.82(5) 


1.84(2) 


1.604(12) 


1.596(11) 


F+- 


2.6(9) 


1.94(6) 


1.950(11) 


1.72(2) 


1.721(8) 



TABLE IX: Comparison of the m3 iPr i or with the results from the full 3 exponential fit and the 
Bayesian fit for "DF3sum6" configurations. Again, the prior 7773 values are in fairly good agreement 
with the ttt-3 results from the full 3 exponential fits. Fixing 777.3—7712 does not change the first excited 
state 7772- See Table IVTHl for notation. 



Lattice 


Direct 


Indirect 


Lattice 


Direct 


Indirect 


Lattice 


Direct 


Indirect 


DF3plain 


0.26(4) 


0.41(12) 


DF4plain 


0.69(4) 


0.75(14) 


DF5plain 


0.50(6) 


0.7(2) 


DF3sum6 


0.15(4) 


0.32(9) 


DF4sum6 


0.07(5) 


0.27(12) 


DF5sum6 


0.13(6) 


0.40(16) 


DF3hyp 


0.00(4) 


0.00(11) 


DF4hyp 


0.03(5) 


-0.06(17) 


DF5hyp 


0.07(7) 


0.0(2) 



TABLE X: P-wave spin-orbit splitting vqAE = ro(7n,ip + — tt7ip_) for the different lattices. To get 
AE in GeV requires a factor 0.38(2). 

to the x 2 . This is not a hard constraint, unlike fixing 777.3 to a given value would be, but rather 
constrains the parameter to a given range softly. The m 3j prior and Am 3j prioT are determined 
beforehand by estimating the difference A77732 = — mi from full 3 exponentials fits. 
This mass difference seems to be almost constant for states that have L = 1 or higher (see 
Fig. |24"1). Therefore we use the P_ state to set the A77732 for D-wave and F-wave states. The 
wis, prior for D-wave and F-wave states is then calculated by adding A77732 to the 7772 from 
the full 3 exponential fit for the state in question (see Tables IVIIIl IIXI) . The prior 7773 values 
are in fairly good agreement with the 7773 results from the full 3 exponential fits, and fixing 
m 3 ~ m 2 does not change the first excited state 777 2 . The P + "Bayes" fit are used to check 
that the Bayesian ideas work well and does not restrict the analysis too much. 

C. Spin-orbit splitting 

One interesting point to note here is that the spin-orbit splitting of the P-wave states is 
small, almost zero, for the preferred "hyp" smearing. We extract this energy difference of 
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FIG. 25: (Color online) The Spin-Orbit splittings of P-wave states for the "DF3" lattice. 
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FIG. 26: (Color online) The Spin-Orbit splittings of P-wave states for the "DF4" lattice. 



Lattice 


Direct 


Indirect 


Lattice 


Direct 


Indirect 


Lattice 


Direct 


Indirect 


DF3plain 


0.13(9) 


0.1(3) 


DF4plain 


0.64(9) 


0.6(2) 


DF5plain 


0.20(10) 


0.2(3) 


DF3sum6 


0.13(3) 


0.13(11) 


DF4sum6 


0.17(6) 


0.2(4) 


DF5sum6 


0.43(5) 


0.5(4) 


DF3hyp 


0.28(5) 


0.27(14) 


DF4hyp 


-0.18(7) 


0.2(4) 


DF5hyp 


0.12(6) 


0.1(5) 



TABLE XI: D-wave spin-orbit splitting tqAE = ro(miD + — wiid_) f° r the different lattices, 
get AE in GeV requires a factor 0.38(2). 
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FIG. 27: (Color online) The Spin-Orbit splittings of P-wave states for the "DF5" lattice. 
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FIG. 28: (Color online) The Spin-Orbit splittings of D-wave states for the "DF3" lattice. 

the 1P_|_ and 1P_ states in two different ways: 

1. Indirectly by simply calculating the difference using the energies given by the fits in 
Eq. El when the P + and P_ data are fitted separately. 

2. Combining the P + and P_ data and fitting the ratio C2(P+)/C2(P_), which enables 
us to go directly for the spin-orbit splitting, mip + — mrp_. 
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FIG. 29: (Color online) The Spin-Orbit splittings of D-wave states for the "DF4" lattice. 
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FIG. 30: (Color online) The Spin-Orbit splittings of D-wave states for the "DF5" lattice. 

In the latter case, the expression (for a given fuzzing) is 
C 2 (P + 



C 2 (P. 



Ae 



-Ami T 



l + bte- Am ^ T + bt e" 



1 + b~ e- 



-Atn; T 



+ h e~ 



(8) 



27 



where 



Ami 


= ™ip+ 


- m 1P _ , 




= ^2P+ 


- ^ip+ , 


/\ TTi 


— 777 on 


'"■IP- ! 


Am^ 


= ™3P+ 


— m X p + and 


AmJ 


= ^3P_ 


- rn 1P _. 



We get the best results by fitting Ami, Am^, Am^ and the coefficients A, and 6J, but 
fixing the remaining mass differences and 6's from the individual two-point correlator fits 
(equation [2] and Table [IT]) • Thus 

= a 3 , /l (P + K /2 (P + ) = a 3Jl (P-)a 3j2 (P-) 

3 a lj/l (P + )a lj/2 (P + ) 3 a lth {P.)a ltfa (P.) 1 J 

for given values of fuzzing indices fx, fi- The D-wave spin-orbit splitting is also extracted 
in a similar manner. The results of the fits are given in Tables |Xl |XI] and in Figs. [231 - [5U1 
In all cases the errors on the direct estimates are much smaller than those on the indirect 
ones. Also in most cases the direct and indirect estimates are consistent with each other - 
the only exception being the P-wave "sum6" estimates. There the direct value is somewhat 
lower than the indirect estimate. In fact this difference brings the "sum6" direct estimate 
closer to the "hyp" value, and lends support to the preferred "hyp" estimate, which in 
all three cases gives a small P-wave spin-orbit splitting (SOS), consistent with zero, for the 
"hyp" configurations. The D-wave spin-orbit splitting (SOS) results are more varied, but the 
"DF3hyp" lattice suggests clearly a positive, non-zero D-wave SOS. However, the "DF4hyp" 
and "DF5hyp" estimates are considerably smaller, becoming negative for the "DF4hyp" . At 
present it is not clear whether this is a lattice artefact due to, say, not being in the continuum 
limit, or that indeed the D-wave results are more dependent on m q than in the P-wave case. 



VI. A MODEL BASED ON THE DIRAC EQUATION 

Since the mass of the heavy quark is infinite, we have for a potential description essentially 
a one-body problem. Therefore, a simple model based on the Dirac equation is used to try 
to describe the lattice data. The potential in the Dirac equation has the usual linearly rising 
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scalar part, 5 sc r, but in addition a vector part b vec r is added. The one gluon exchange term, 
aoGE • Voce, where 

v - 4a * (r) nm 

^OGE — , l-WJ 

3 r 

with the running coupling constant a s (r) given by 

. . 2 r„ sin(A;r) 
a s (r) = - / dk— K -^a s (k 2 ) (11) 
Jo k 

and 

2 \ 1 

° s{k } = "27ln[(^ + 4m2)/(A^ CD )]- (12) 
Here, guided by fits to various meson masses using the Blankenbecler-Sugar equation, we 
fix Aqcd = 260 MeV and the dynamical gluon mass m g = 290 MeV (see 19j for details). 
The potential also has a scalar term mwL{L + 1), which seems to be needed to increase the 
energy of D-wave states. This type of term arises in flux tube models, where a flux tube's 
rotational energy is proportional to L(L + 1) (like in Isgur-Paton flux tube model, 2o|). 

The lines in the energy spectrum plot (Fig. l3~Ti) show three Dirac model fits from Table IXTTl 
with m = 560 MeV (the constituent quark mass, from 3]) and ooge = 1.00. Attempts to 
also vary cioge easily lead to instabilities. The solid line, labelled "fit 1", is a fit to three 
"DF3hyp" energy differences: 1P_ and 1D_ with respect to the ground state, and the P- 
wave spin-orbit splitting (direct estimate) SOS(IP) [i.e. E(1P + )— E(1P_)]. The fit to these 
energies is acceptable with total \ 2 — 1-68, but as soon as a fourth state [e.g. SOS(D)] is 
added a good \ 2 can no longer be achieved. The dashed line, "fit 2", shows an attempt to 
fit "DF3hyp" 1P_, SOS(IP), 1D_ and SOS(ID). The X 2 is not good, and letting a OG E vary 
does not help: that only leads to unphysical values for the parameters. Using a different 
constituent quark mass, say m = 490 MeV from [21], gives basically the same fits (the 
changes are minimal). "Fit 3" is a fit to "DF3sum6" 1P_, SOS(IP), 1D_ and SOS(ID), 
and is shown in the figure for comparison. The fits to "DF3sum6" energies are also shown 
in Table Kill In Fig. [32] the same Dirac model fits are shown for the excited states. Here 
it can be seen that the fit is about 500 MeV lower than the lattice results, and the shift 
seems to be constant for both lattices ("DF3sum6" and "DF3hyp") for all states, except the 
2S. There is no obvious reason why the Dirac model should underestimate the first radial 
excitations by a constant amount, but a term of the form 0.5(n — 1) GeV could be included 
in the model to improve the fit to excited states and be interpreted as a flux tube effect in 
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Case 


6 SC [GeV/fm] 


b vcc [GeV/fm] 




total x 2 


sum6 P_ 


SOS(IP) 




1.168(2) 


0.0 


0.0 


7.8 


sum6 P_ 


SOS(IP) 




0.86(13) 


0.57(18) 


0.0 


0.0 


sum6 P_ 


SOS(IP), D_ 




0.372(14) 


0.0 


0.0916(14) 


4.5 


sum6 P_ 


SOS(IP), D_ 


SOS(ID) 


0.265(4) 


0.571(13) 


0.0696(11) 


2.1 


hyp P_, 


SOS(IP) 




1.294(10) 


0.0 


0.0 


0.54 


hyp P_, 


SOS(IP) 




1.40(18) 


-0.2(2) 


0.0 


0.00 


hyp P_, 


SOS(IP), D_ 




0.763(9) 


0.0 


0.0554(9) 


1.83 


hyp P_, 


SOS(IP), D_, 


SOS(ID) 


0.48(5) 


0.26(2) 


0.066(5) 


34.3 



TABLE XII: Dirac model fits for "DF3" . Here ooge = 1 an d constituent quark mass m = 560 MeV. 
Fits are attempted for "DF3sum6" and "DF3hyp". A "perfect" fit (2 fit parameters, 2 data points) 
can be found for P_ and the P-wave spin-orbit splitting, if both scalar and vector linear potentials 
are used. However, all P- and D-wave data [P_, SOS(IP), D_, SOS(ID)] can not be fitted using 
the two linear rising potentials and adding a scalar term muL(L + 1) still does not give a good % 2 - 

the same philosophy as the uL(L + 1) term. However, as the fit to the ground state energies 
is poor, this improvement is not pursued. 



VII. CONCLUSIONS 



With the "DF3hyp" lattice, our predictions for the 1 + and 2 + P-wave state masses 
agree very well with the experimental results. We also predict that the masses of the 
two lower P-wave states (0 + and 1 + ) should lie only a few MeV below the BK and 
B*K thresholds respectively. 

Also with the "DF3hyp" lattice, the P-wave spin-orbit splitting is small (essentially 
zero), but the D-wave spin-orbit splitting is clearly non-zero and positive. In con- 
trast, another lattice group finds the P-wave spin-orbit splitting to be positive (about 



35 MeV) and the D-wave SOS to be slightly negative (see [151]). i.e. they seem to 
observe the famous inversion 21] . However, the recent European Twisted Mass Coll- 
aboration results find the P-wave SOS to be negative and the D-wave SOS to be small 
l| . One clearly needs to go to the continuum limit before any definite conclusions can 
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FIG. 31: (Color online) Energy spectrum of the heavy-light meson and three Dirac model fits. "Fit 
1" is a fit to "DF3hyp" 1P_, SOS(IP) and 1D_, whereas "fit 2" is an attempt to fit "DF3hyp" 
1P_, SOS(IP), 1D_ and SOS(ID) (see Table [XTH). "Fit 3" is a fit to "DF3sum6" 1P_, SOS(IP), 
1D_ and SOS(ID), and is shown here for comparison. 
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FIG. 32: (Color online) Energies of the first radial excitations of the heavy-light meson and the 
same Dirac model fits shown in Fig. [3TJ 
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be made. 

In [ijj] Woo Lee and Lee suggest that the absence of spin-orbit inversions can be ex- 
plained by chiral radiative corrections in the potential model. Small spin-orbit split- 
tings throughout the meson spectrum could be explained by a relativistic symmetry in 



the Dirac Hamiltonian discussed in 



221 ]. This would indicate that the scalar potential 



is (at least approximately) equal to the vector potential. 

The one-body Dirac equation model with one-gluon exchange, vector and scalar linear 
potentials and a scalar term muL(L + 1) (like a flux tube rotational energy) is not 
good enough to describe the entire lattice energy spectrum. Therefore, one should be 
very careful in using such simple potentials to describe the interaction between quarks. 
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APPENDIX A: SOME CHECKS 



To check that the different results for different smearings (specifically for the spin-orbit 
splittings) are a real effect and not due to some biasing element in the final analysis, we plot 
the basic signals in Figs. ESI for the P-wave. The signals clearly are different for "DF3sum6" 
and "DF3hyp" , supporting the results of the complete analysis in Fig. [23 
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